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Abstract 

We exhibit the limit shape of random Young diagrams having a distribution pro- 
portional to the exponential of their area (grand-canonical ensemble), and confined in 
a rectangular box. The Ornstein-Uhlenbeck bridge arises from the fluctuations around 
the limit shape. The fluctuations for the unconfined case lead to a two-sided stationary 
Ornstein-Uhlenbeck process. 
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Introduction 

A partition of an integer n is a finite non-increasing sequence of integers 

7T = (7Ti,7r 2 , . . . ,7T fc ), 

with 7Ti + • • • + 7T/% = n. A convenient graphical representation in (Z + ) 2 of the partition 7r 
is its Young diagram, also denoted by ir. The Young diagram consists of a stack of 7Tj unit 
squares on the i-th column. See Figure [TJ The area \tt\ of the Young diagram ir is the total 
number of squares, which is equal to n. 

If the number of summands k is less than some integer a and all the summands ttj are 
less than b, then the Young diagram stays in the box [0, a] x [0, b]. 

In this paper we study the limit shape and fluctuations of large random Young diagrams 
assigned to stay inside a (large) rectangular box, when the probability measure is propor- 
tional to q , l 7r l, and q is a suitable parameter. It is easy to see that if one doesn't want the 
system to degenerate in the limit, one has to make q depend on the size of the box. Namely 
1 — q has to be taken of the order of the inverse of the side of the box. 
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Figure 1: The Young diagram for the partition ir = (8, 7, 7, 5, 4, 3, 3, 2, 2, 1) of n = 42. 



The study of the combinatorics of partitions of integers goes back to Hardy and Ra- 
manujan and then Erdos in the 40's. The statistical physics point of view was introduced by 
Vershik [9] in his study of the typical shape of the partition of a large integer and obtained 
what we will call Vershik's curve. Recently, Funaki and Sasada obtained Vershik's result, 
in [5], as a by-product of an hydrodynamic result for the corresponding particle system 
model. 

Our study lies at the intersection of classical topics of probability theory, combinatorics, 
and statistical physics. 

Along the paper, we make an extensive use of the classical Gauss polynomials, well- 
known in combinatorics as the generating functions for the number of Young diagrams 
with given area (see e.g. [I]). In order to get to the asymptotic regime, we are led to state 
a g-analogue for Stirling's formula. 

For all values of the parameters of the problem (limiting aspect ratio of the box, and 
c = limn(l — q)), the limit shape obtained turns out to be a restriction of Vershik's curve 
as it was recently noticed by Petrov in [6]. 

The main part of the paper deals with the fluctuations around this limit shape and 
requires a fine understanding of the boundary of the random diagram. From a probabilistic 
point of view, the boundary is described by a Markov which can be considered as the q- 
analogue of the bridge of the simple random walk. In the classical case (c = 0), the properly 
rescaled interface converges to the Brownian bridge. In the general case, the limiting process 
is an Ornstein-Uhlenbeck bridge. Let us mention that, in the framework of unlimited 
partitions, and under some specifications on the summands, Vershik and Yakubovich [8| II 1| 
already observed Gaussian fluctuations. For strict partitions, see [10J. The fluctuations for 
the unconfined case lead to a two-sided stationary Ornstein-Uhlenbeck process. However, 
in contrast with the limit shape problem, there is no direct argument to deduce fluctuations 
for partitions in a box from the unlimited case. 

The paper is organized as follows. After a presentation of the combinatorics of the 
problem (Section 2), we perform its asymptotic analysis through a q- version of Stirling's 
formula (Section 3). Section 4 is devoted to the derivation of the limit shape phenomenon. 
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In Section 5, we study the fluctuations around the limit shape: we first compute the limiting 
2-correlation function and deduce the convergence of finite-dimensional marginals. This 
section ends with the delicate proof of the tightness of the fluctuations. 



1 Presentation of the model 

We study the asymptotic distribution of random Young diagrams fitting in a large rectan- 
gular box with dimensions a x b. Given a real number q > 0, we assign to each diagram it 
the probability 

where \ir\ is the number of boxes of tt and Z a ^(q) is the partition function, the sum of all 
qfr\ _ For the sake of clarity, the probability measure P^ b will be simply denoted, by P. 

We fix a parameter p 6 (0, 1), and choose the dimensions of the box for each n to be 
a n x b n , where (a n ) and (b n ) are sequences of positive integers satisfying: 

a n + b n = 2n, lim = - . 

n->-oo n 1 — p 

We are interested in the limiting behavior when n goes to oo. To obtain a non-degenerate 
limit, q must go to 1 as n goes to infinity. We fix a real parameter c and pose q = e _ «. 

Note that c = corresponds to a uniform probability measure. We will assume that 
c 0. The results for c = can be obtained by taking limits. The physical meaning of the 
parameter c is that of a pressure, since it is the variable conjugated to the (two-dimensional) 
volume. 



2 Combinatorics of the partitions 

Let us start with a fixed box with dimensions a x b. The partition function Z a ^(q) is 
expressed in terms of Gaussian polynomials or (/-binomial coefficients, where integers j are 
replaced by their (/-analogues (j) q = ^E^". 

Definition 1. Let q > 0, 



n n 

^=n^=n T 



q> 



i=i" Q 

n\ n\ q 1 - q n ~ j 



mj q {n — m)\ q m\ q j-J^ 1 — q m i 
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Lemma 1. For all (a, b) £ (N*) 2 , and all q > 0, the partition function Z a ^(q) is equal to 

^K°:V(°r), 

Proof. ( a + 6 ) and Z a<b (q) both follow the recursion relation 

z a ,b{q) = z a ,b-i{q) + q b z a -i,b(q), 

where the first term corresponds to those diagrams with all parts strictly smaller than b, 
and the second term to the diagrams with at least one part of size b. □ 

We use the following coordinates. The bounding rectangle is the rectangle in the plane 
with corners at the points (0, 0) and (a + b, b — a) and sides with slopes ±1. The boundary of 
a diagram is encoded as a lattice path {Xk) 0<k<n from the origin to the point (a + b, b — a) 
and such that X^+i = ± 1 for all k. 




Figure 2: A partition as a lattice path. 

Computing the probability P {X^ = t) is straightforward. A lattice path passing through 
(k,£) is composed of a path from (0,0) to (k,£) and a path from (k,£) to (a + b, b — a). 
Since X k has the same parity as k, the exact formulas also depend on this parity. 

Proposition 1. The 1- dimensional marginal of X under F is given by 

q(k+i)(a~k+i) 

PoftP^fc = 2i) = — —- — Z k _i^ +i (q)Z a _ k+i £_ k _i{q) (1) 

^a,b(q) 
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and 

q(k+i+l)(a-k+i) 

^ q a h{X2k+i = 2i + 1) = — -r-: Z k _i tk+i+ i(q)Z a _ k+ i ib _ k _i_i(q) (2) 

z a ,b{q) 

We focus our attention to the behaviour of the process at even times, which is sufficient 
to study its scaling limit, since the process has bounded steps (+1/— 1). 




k — i 



Figure 3: Illustration of Equation ([I]). 

We now mention a unimodality result for the distribution above which turns out to be 
useful later. 

Lemma 2. The function I \— > ¥ q ab (X2 k = 21) is unimodal: there exists an integer L^(k) = 
l) a b (k) such that 

P£ 6 (*2* = 2(* + l)) 



< i 



t > i){k). 



(3) 
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Proof. Writing the ratio 

K,b( X 2k = 2(1 + 1)) = 2 , fc _^_i ifc+ ^ + i(g)Z a _ w+ i i b_ fc _^_i(g) ^ (a _ fc+ ^ +1)(w+11 _ (a _ w)(fc+< ; ) 

This ratio is smaller than 1 if and only if 

(1 _ q k-i ){ l _ q b-k-t )q a+2t+l < (1 _ _ 

^ ( g a+6+l _ 1) < _ l) g 2£ + (^ q k+l( q a _ x) + _ ^ ^ 

Both terms on the right hand side are increasing continuous functions of I, proving the 
existence of an integer L$(k) with the asserted property. □ 

Proposition 2. The 2- dimensional marginal of (^2fe)o<fc<n is given by 
Fl b (X 2k = 2i,X 2l = 2j) = 

Zk+i,k-i{q)Zj+l-i-k,l+i-k-j{q)Z a -l + j,b-l-j{q) ^g~(k-i))(k+i)+(l+j-k-i)(a^l+j) 

z a ,b(q) 

Now we have exact expressions for probabilistic quantities, we let the dimensions of 
the box depend on n and investigate the limiting behaviour of the random Young diagram 
under P-? , in the regime 

a n + K = 2n — > +oo, — — > p G (0, 1), — n log q — > c G M. 

2n 

3 Asymptotics of q-factorials 

In order to study the limit of the process as the size of the box goes to infinity, we first 
need the asymptotic behavior of the ^-factorial. 

Proposition 3. Let e > 0. For I > (1 — q)~ 1 e, 

(£\) q = ^2^exp QT log^d^ (l + 0(l-g)) 

beginequation* Consider log£! g == X^=il°g^<?- By the Euler-Maclaurin formula, we 
have 

log£! g = ^ log(gdt + i(log^) + y Bi^—g-i-y Iog(t,)dt. 
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where B\ (t) = t — [t\ — \ is the first periodic Bernoulli polynomial. 

The integral \ogt q dt = 1 + 0(1 — q). For the second integral in ([3]), we add and 



subtract ~ to d tq in order to get 



(4) 



The function /(i) = d ^ f 9 — f is continuous on and its derivative is bounded by an 
absolute constant M times (1 — q) 2 , i.e. \f'(t)\ < M(l — g) 2 for all t G R+. On the interval 
[A; — l,k], the function B\{t) is given by = i — k + ^ and has the antiderivative 

Bn,k(t) = \t 2 ~ (k — \)t + fc( ' fc 2 ~ 1 ^ with B n j t (k — 1) = B n ^(k) = 0. By partial integration 
we get 



/ Bx(t)f(t)dt < MV / |Bn,j|,(t)|dt 



M£(l - g) 5 



12 

fe=l u K 

proving that this term is 0(1 — q) asn-f oo. The first term of Q is found to be 
J B l {t)\dt = £-l-{£ + \)\og£ + \og£\ = -1 + \ log 2tt + 0(±) 



(5) 



by the classical Stirling approximation £\ = \Jlixl (^) (l + 0(|)). Combining ([3]), (jl]) and 
([5]), we get the result. 

Using the notations of Section Q] we restate the q-Stirling's formula in a convenient form 
for this particular context. 

Corollary 1 (q-Stirling's Formula). Let c £ and fix e > In the limit when n goes to 

c 

infinity, with q = en, the following asymptotics hold for all £ > en 



where 



% = V2^\ ^—±n £ eMnS c (i)){l + 0(f)), 



S c {a) = J log I J dx = aS ac (l) 



In particular, there exist positive constants m and M such that for all n, and all £ 
between 1 and n, 



iV2^n~\ C " - V exp(nS c (|)) < £\ q < mV2^\ 6 " - - n e exp(ra5 c (^)). 
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We turn now to the the asymptotics of the 1- and 2-dimensional marginals of X under 
given by Propositions Q] and [2j For this purpose, define 



f c (x, y) = S c (x + y) - S c (x) - S c {y) and h c (x, y) 



' (e cx - l)(e c y - 1)' 
If x = £ and y = ~, the asymptotics of the g-binomial coefficient is given by 



Z l,k 



As a consequence, 



IS = ( i+ £ k ) = ^M^y)exp(n/ e (x,y))(l + 0(l)). 



Corollary 2. For all e > 0, for all k between en and (1 — e)n, with s = k/n and x = i/n, 
we have 

F(X 2k = 2i) = J ^-H$ (s, x) exp(nF { V c ( S , x)) (l + 0(1)) , 



2im 

and there exists a constant M > such that for all n, i and k, 

P(X 2fc = 2i) < My^i7W( S ,x)exp(nF p ( 1 c )( S , a ;)), 

where 

x) = -c(2p - s + x)(s + x) + / c (s - x, s + x) + / c (2p - s + x, 2 - 2p - s - x) - / c (2p, 2 - 2p), 
/i c (s + x, s — x)h c (2p — s + x, 2 — 2p — s — x) 



H$(s,x) 



h c (2p, 2 - 2p) 



In the first equality, when replacing |? by p in the indices of .H^ 1 ) , we make an error of 
order O(-), which is absorbed in the factor (l + 0(— )). Making the same substitution in 
the indices of F^ would change the multiplicative constant in the asymptotics. 

Corollary 3. For all e > 0, for all k < I between en and (1 — e)n then, with s = k/n, 
t = l/n and x = i/n, y = j/n 

P(X 2k = 2i,X 2l = 2j) = ^-H^(s,t,x,y)exp(nF^ c (s, t, x, y))(l + 0(1)), 

Inn ' 2n > c " 

where 

F® (s, t, x, y) = -c((2p -s + x)(s + x) + ( y t-s + y-x)(2p-t + y)) 
+ f c (s - x, s + x) + f c (t - s + y-x,t-s — y + x) 
+ f c (2p - t + y,2-2p-t-y)- f c (2p, 2 - 2p), 
(2) h c (s + x,s - x)h c (t - s + y - x,t - s - y + x)h c (2p - t + y,2 - 2p - t - y) 

H pM^y) = h c (2p,2-2 P ) • 
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We give now a slight refinement of Corollary [3] that will be useful in the fluctuations. 



Corollary 4. For all e > 0, for all k < I between en and (1 — e)n then, with s = k/n, 
t = l/n and x = i/n, y = j/n, for all s' = s + o(l), t' = t + o(l), x' = x + o(l), y' = y + o(l): 



2LX- 



21 



2j) 



1 



Hf}{s\ t',x', y') exp(nF® (s, t, x, y))(l + o(l)). 



2-Kn 



4 Limit Shape 

We associate to the lattice path (X^) the continuous piecewise linear function. X : s t— > X s 
defined on [0, 2n], which coincides with when s = k. The graph of the function X is the 
boundary of the random Young diagram we consider. 
Let L PjC be the function on [0, 1] defined by 

p i , . 1, sinh(ct) + e^- 1 ) sinh(c(l -t)) 

Vt € 0, 1 , L p>c (t) = 1 - 2p + - log ^ — ^ U. 

c sinh c 



7 ^g 



e ct — e ct + e" 



-c(2-2p-t) 



-c(t-2p) 



-c(2-2p) _ p 2cp 



(6) 



This is the limit shape of the rescaled random process, in the sense made precise in 
Theorem[T]below. First we need to estimate the proximity between the value of the function 
L P)C and the most probable value for denoted by Ln(k) = L an bn e -c/ n (k) in Lemma [2l 



Lemma 3. For all n and allO < k < n, 



±4(*) - < 1+ | - P\ , ™here q = e^. 



Consequently, in the limit when n goes to infinity, and a n /2n goes to p, 



o. 



(7) 



Proof. The ratio of probabilities ([3]) used to define L${k) = Ln(k) can be rewritten in terms 
of the function 



(p,c,t,x) (->■ R pc (t,x) = e" 



,sinh(|(t - a?)) sinh(§(2 -2p-t-x)) 
sinh(§(t + x)) sinh(§(2p - t + x)) 



namely as R±n c (-, -). The function x \- > R pc (t,x) is decreasing, and by definition L„(/c) 
is the smallest integer £ such that Ra n ( — ,—)< 1. 

2n ' 

On the other hand, a computation shows that for each t £ [0, 1], the equation 

#/9,c(*, x) = l 

has the (unique) solution x = L P}C (t). Substituting |^ for p, we conclude that 



^L{(k)-La n (h.) 



2n' 



< I 



(8) 
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for all n and all k < n. 

Differentiating L PjC with respect to p, we find that 



dp 

for all t, and so by the mean value theorem, 



< 1 



2n ' 



< 



Yn~ P 



for all k < n. From (El) 



(9) 
□ 



j) and the triangle inequality, we get the expected result. 

Remark 1. It turns out that log R p c (t, x) coincides with the derivative of Fpl}(s,x) with 

respect to x. As a consequence, L PjC (s) can also be viewed as the argmax of Fp^}(s, •), which 
is nonpositive and vanishes at L P)C (s). See Corollary [6j 

Theorem 1. The boundary of the rescaled random Young diagram converges in probability, 
for the uniform topology, to the curve of t ^ L p ^ c (t). 



Ve > 0, lim roe ° 



(in ,b n 

Vte[o,i] 



sup \^X 2tn 



L p , c (t)\ >e 1=0, 



where 



L p Jt) = Mog 



-ct 



e ct + e 



-c(2-2p-t) 



-c{t-2p) 



e -c(2-2p) _ e 2cp 

Remark that in the case of a square box (p = 1/2), the expression for L PjC boils down 



to: 



1 i cosh(c(t - §)) 
c ^ cosh £ 



Proof. Fix e > 0. For t < | or i > 1 — | , the difference | ^^tn — L/9, c(t) | is always smaller 
than e. We have to control what happens for t E (|, 1 — |). Using the fact that L PiC is 
differentiable, and its derivative with respect to t is bounded by 1, we have that 

1 



\L P ,c^[tn\) -L p>c (t)\ < 



n 



which is smaller than | for n sufficiently large. The same is true of 
Thus, by an |-argument, to control the sup over (|, 1 
happens at points of the form t = —: 



i 



Xo 



|), it is sufficient to control what 



SU P I 2n 

v*e(§.i-f) 



1 X2tn 



L P,c{t)\ > £ 



< Yl P ( X2fc > MLpATi) + ° r x 2k < 2n(L P)C (|) - e)) (10) 



fceNnn(f ,i-§) 
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But it follows from (J7|), that for n sufficiently large, 2n(L(^) + e) > 2Ln(k) for all k € 
N n 1 — |), and thus using the unimodality of the law of X 2k 

F(X 2k > 2n(L AC (|) + ej) = £ F(X 2k = 21) 

«>n(L p , c (|)+ e ) 

<nxF(X 2k = 2[n(L PjC (^) + e)\). 

which by Corollary [2] is exponentially small, uniformly in k, as n goes to infinity. Thus for 
n sufficiently large, the sum on the RHS of ([TO]) can be made smaller than any positive 
number. □ 

Remark 2. From the proof and Lemma [3] we see that the convergence statement in the 
theorem holds uniformly for a family of sequences (a n ) and (b n ), so long as the convergence 
Tjr- — > p is uniform for the family. 

5 Fluctuations 

We now study the fluctuations of the interface around the limit shape. We define for n 6 N* 
a new rescaled process: 

Vie [0,1], x t = x[ n) = (J^x 2nt - L P)C {t)\ . (11) 

We place ourselves in the space D of cad-lag paths on [0, 1] endowed with its usual 
topology. We state now our second main result for the convergence of the fluctuations of 
the interface to the Ornstein Uhlenbeck bridge. See the appendix for the definition and 
some properties of this process. 

Define 



_ 1 y/(e 2c P - 1)(1 - e- 2c d-p)) _ ^/2sinh(cp) sinh(c(l - pfj 



^2 sinh cs + e^P- 1 ) sinh(c(l - s)) e <\~p) s inh(cs) + e c( - p ~^ sinh(c(l - s)) ' 

(12) 

Theorem 2. The sequence (X^ / f(s)) n , converges weakly in D to the Ornstein-Uhlenbeck 
bridge (Xs) se [o,i]; which is the Gaussian process on [0, 1] with covariance 

sinh(cs) sinh(c(l -t)) 
csmh(c) 

for < s < t < 1. 
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5.1 Two-point correlations 

To prove convergence of the two-dimensional marginal to a Gaussian process, we apply 
a saddle-point method, ie. we need to show that the function F^ 2 \ which governs the 
exponential decay of 2-dimensional marginals (see Corollary [3|) , has a critical point 'on the 
limit shape', and that it takes its maximal value of at this point. We will prove that this 
is indeed the case, but first we need a lemma describing the limit shape in a subrectangle. 

Lemma 4. The limit shape L satisfies the relations 



and 



for < s < t < 1, with 



-^- s {L p , c {t) ~ L p , c (s)) = V,c(fE§) (13) 
h p>c (s) = Vv (f) (14) 



and 



, 2p- s + L PjC (s) , 
P= 2(1-,) ■ c=c(l-s) 

= t - L PtC (t) „ = 

p 2t 

Proof. These relations can be checked analytically. We sketch now a less computational 
argument for Equation (|14p . The same strategy applies for Equation (I13p . Fix t £ (0, 1), 
and take a sequence of boxes with sidelengths [0, a n ] x [0, b n ] (with a n +b n = 2n, a n /2n — > p). 
Theorem Q] states that, when n goes to infinity, the probability under c b " that X n t/n 
converges to L PiC (t) goes to 1. As a consequence, the restriction of the limit shape L P)C to 
the, rescaled by n, lower left subbox [0, LJ^dH j x [q, * +l p- c W ] i s the rescaled limit shape 
for boxes with ratio t L ^v^^ — = 1 L £f^ ■ The rescaling factor is the perimeter of 

/ -dL 

the subbox t. The value of the parameter q remains the same. Since q = e c ' n = e «* , we 
get c" = tc. □ 

This lemma basically describes the limit shape of the process restricted to a subrectangle 
defined by a point (s, L PtC (s)) and either the right or left corner of the original rectangle. 

Lemma 5. For all < s < t < 1, the functions and F^ satisfy the relations 

F$(s,t,x,y) = F$(s,x) + (1 - s)F^ c ,(t',y') (15) 

and 

F$(s,t,x,y) = F$(t,y)+tF { V( S ",x") (16) 
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with 



and 



, 2p-s + x t-s , y-x , 

P = — ; ^! t = , y = , c = c(l — s) 

H 2(1 -a) l-s 1-s V ' 



n t y „ s „ x 



x" = -, c" = tc. 



2t. r t 

Proof. The relation (|15p follows from the identity 

F(X 2k = 2i,X 2e = 2j) = F(X 2k = 2i)F(X 2e = 2j\X 2k = 2i). 

allied to scaling arguments similar to those used in the proof of the previous lemma. Equa- 
tion (fT6|) follows from the conditioning on X 2 i instead of X 2k . □ 



Lemma 6. The partial derivatives of with respect to p, s and x all vanish at x = 
Lp, c {s). 

Corollary 5. The partial derivatives of with respect to p, s, t, x and y all vanish at 
x = L P)C {s) and y = L P)C (t). 

aF( 2 ) i aF< 2 > 



Proof. For the partial derivatives with respect to p, use (|16p . For ds and Qx , use ([16 
together with ([II]). Similarly, for and ^g^-, use (JEJ) together with (JTHJ)- C 

Corollary 6. F$ (s, L p>c (s)) = and F$(s,t,L PiC (s),L PiC (t)) = 0. 
Proof. By the chain rule and Lemma [6] we get 



ds 



(F$(s,L p , c ( S )))=0. 



Taking the limit s — )■ yields zero, proving the claim. For the second identity, use the first 
one together with (fTUj) and ([Tl"|) . □ 

Proposition 4. Assume ^ = p + o(-^) . Let s < t in [0, 1]. T/ie joini lam o/ (xi n) , X t (n) ) 
converges to the law of the centered 2- dimensional Gaussian vector 

(X s ,X t ) = (f(s)Y s ,f(t)Y t ), 

where f(s) is defined in as in Equation (|12p . and (5^)fg[o,i] is an Ornstein-Uhlenbeck bridge 
on the interval [0, 1] with parameter c (see Appendix). 

Remark 3. In the case of a square (p = i), the expression of / simplifies drastically to 

1 



v^cosh (c(s - §)) ' 
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Proof. Let ai, 61,02, 62 be four real numbers such that a\ < b\ and a 2 < 62- In the following, 
we write L for L PjC . Clearly, X^ G [ai, 61] if and only if X 2ns is in the interval 

2nL(s) + [2ai-v/n,26i-v/n] . 

Since |X2 ras — -X"2|_nsJ | < 2 and |X2 n j — -X~ 2 |_ nt j | < 2, we just need to compute the limit of 

W((X 2[nsi ,X 2[nti ) 6 Ai x A 2 ) ■ 

as n goes to infinity. 

Set k n = [ns\ and i n = [nt\ . One has: 

P ((X 2LnsJ , X 2[nt} ) £ A 1 xA 2 )= P ^ = 2i, X 2 , n = 2j) 

(i,j)e| J 4ixiA 2 

(i,j)e±Aix±A 2 

Here, i 7 ^ 2 ) < everywhere and the number of terms is O(n), so the error terms o(l) 
can be replaced with an o(l) term outside the sum. Thus, we are left with a Riemann sum 
for the double integral 



2im 



H$(s,t,L(s),L(t)) / / exp 



«41(^T>^)ld^ (17) 



in which we make the substitution 

u = x\/n + nL(s), v = y\fn + nL(t) 

to get 

±H$(s,t,L(s),L(t)) [[^V [nFg c ^/ t ,^ + L(s),^ + L(t))} dydx 

We make a second degree Taylor expansion of the function Fa2 (— j t) a ^ the point 
(L(s), L(i)). Let z„ be the point L(s),L(t)) in M 4 . 

nF¥>Jz n ) + ^-(z n )xV^+^-(z n )y^ (i 8 ) 



+ 



<9x <9y 



2 <9x 2 



^(z„)x 2 + 



2 9y 2 (^ 2 + ^(^ + 0(n-,) 
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We must find the limit of this expression as n — > oo. The first term on the right has 
limit zero, which can be seen by Taylor expansion around (s, t, L(s), L(t)), using Corollaries 
E and E 

Since is an analytic function, we have in general 

F m , (s + o(n-^),t + o(n-^),x,y) =F$(s,t,x,y)+o(n-^) (19) 

and similarly for all partial derivatives of F^ 2 \ Since = p + o(n~2) ; ^ = g + o(n~z) 
and ^ = i + o(n~2) 5 this shows that the second and third terms on the right in (j 1 8|) tend 
to zero as n — > oo. Thus, if z = (s, t, L(s), L(t)), 



ld 2 F t 



217(2) 



(20) 



as n —7- 00. From (|19p it also follows that the sequence F^_ (%■)%>•>') i s equicontinuous, 

2n ' 

hence by Ascoli's theorem that the sequence converges uniformly, hence that the double 
integral (|17p converges to 



±H$(s,t,L(s),L(t)) 



b r d 



cxp 



1 r)2 77(2) 

- p '° (z)x 2 + 

2 &r 2 {Z)X + 



ld 2 F t 



(2) 



277(2) 



p.c 



2 dy 2 



(z)y 2 + 



d 2 F 



p.c 



dxdy 



{z)xy 



dy dx 



To find these double derivatives, we use (115|) and (fTBT) . Hence we need 



2 77(1) 



d F, 



p.c 



2c\„2cs 



{s,L{s)) = c (l-e- zc )e 



dx 2 

and exploit the fact that 
dx 2 



(1 _ e -c(2-2p) _ e -2cs + e -c(2 S -2p))2 

(1 - e 2c ^)(l - e- 2cs )(l - e - c (2-2p))(l - e - c ( 2 -2s)) 



(s,t,L(s),L(t)) = - 



1 fl2 P (l) 

(s ,L p // c »(s )) 



a2 77(2) 

* ' AC -( S ,t,L( S ),L(i)) 



i <9x 2 
1 <9 2 F (1) 



<9y 2 ' ' ' 1 — s dx 2 

(2) 

This way, we find the double derivatives of FpJ{s,t,-,-), evaluated at the critical point 
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(x,y) = (L(s),L(t)), to be 

d 2 F (2) c(1 _ e 2ct )(1 _ e 2cs + e -c(2-2p-2s) _ e 2c P) ; 



-(s,t,L(s),L(t)) 



dx 2 v '' w ' wy (l-e 2cs )(e 2rf -e 2cs )(l-e 2 ^)(l-e- c ( 2 -^)) 

— csinh ct 
sinh cs sinh c(i — s)f(s) 2 ' 

a 2 F p (2 c } . + r , w ... c(e 2cs - e c )(l - e 2ct + e -( 2 - 2 P- 2 *) - e 2 ^) 2 
-{s,t, L{s),L{t)) 



dy 2 v '' w ' wy (e 2cs -e 2rf )(e c -e 2c *)(l-e 2 ^)(l-e- c ( 2 - 2 ^)) 

— csinh c(l — s) 
sinh c(l — t) sinh c(t — s)f(t) 2 

8 2 F$ c(l — e 2cs + e- c ( 2 ~ 2 ^ 2s ) - e 2c Q(l - e 2rf + e ~ c ( 2 ~ 2 ^ 2t ) - e 2c P) 

-^-( S ,t,L(s),L(t)) - (e 2 CS _ e 2ct)( 1 _ e 2 C p) (1 _ e -c(2-2p)) 

c 



sinhc(t-s)/( S )/(t) 

The last one, the mixed derivative, was calculated from scratch, i.e. by calculating 

Q2p(2) (2d _ 2cs\ 

dFp ' C -(s,t,x,y)- 0(6 6 } 



dxdy K ' ' ^ ( e c(t+») _ e c(s+x))( e c(t-y) _ e c(s-x)) 

and evaluating at the critical point (x, y) = (L(s), L(t)). Thus, the Hessian of F^f)(s, t, •, •) 
at the critical point is 



p,c sinh cs sinh c(l — t) sinh c(t — s) 

j^y \ / sinhci sinhc(l — t) — sinh cs sinh c(l — i)\ [ 
Jpy i y— sinh cs sinh c(l — t) sinh cs sinh c(l — s) 



The covariance matrix of the limiting Gaussian distribution is the negative of the inverse 
of H(F^), which we compute to be 

_ 1 ff(s) \ /sinh cs sinh c(l — s) sinh cs sinh c(l — t)\ ff{s) 
csinhc \ fit)) \sm~h cs sinh c(l — t) sinh ct sinh c(l — i)/ \ /(*), 

(21) 

The matrix in the middle together with the factor — ir— is the covariance matrix of the 
Ornstein-Uhlenbeck bridge on the interval [0, 1] with parameter c (see Appendix). Further 
computations reveal that 

; 1 = HW(s,t,L(s),L(t)) 
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since both sides equal 



sinh c 



f(s)f(t) y sinh cs sinh c(t — s) sinh c(l — i) 
This completes the proof. □ 

We explain now that for m points < t\ < ti < ■ ■ ■ < t m < 1, the limit of the 
corresponding m-dimensional marginal [X^ ,X^\ . . . , 3tj^ ) is Gaussian, with covariance 
matrix E = (cTj ,) defined by 

Eij = — — — f(ti)f(tj) sinh cti sinh c(l - i,-) 
csmhc 

for all i < j. 

A first approach would be to repeat similar computations as for the 2-dimensional case, 
and observe that the matrix a is the inverse of the m-dimensional Hessian of {t\ , . . . , t m ) . 

Another way takes advantage of the Markov property satisfied by the process (X n ). 
Denote by P^{x,y) the transition kernel of this Markov chain. The law of the m-tuple 

ft 1 *-'-^') - g-en by 

Pot (O^iKlO^) • • • P^ 1>tm (* m -i,x m ). 

Proposition |4] implies the convergence of each ^/nP^^. + _ ^ (x n ,y n ) to the density of the 
kernel of the Ornstein-Uhlenbeck bridge as n — > oo, and (x n ), (y n ) converge. The conver- 
gence in law follows from Lebesgue dominated convergence theorem. Domination is ensured 
by Lemma [8] below and unimodality stated in Lemma [2j 



5.2 Tightness and proof of Theorem [2] 

We already proved the convergence of the finite-dimensional distributions. We need to show 
that the sequence of the distribution is tight. For this we use Theorem 13.5 p. 142 of [2j. 
The criterion is checked in Lemma [TU1 

Before entering the proof, we need a geometric definition. 

Definition 2. Let ABCD be a rectangle with sides having slopes ±1 in the (s,x) coordi- 
nates. For e G (0, 1), the e-parallelogram of the rectangle ABCD is the unique parallelo- 
gram with diagonal AC and sides with slopes 1 — e and — 1 + e in the (s,x) coordinates. 
The e-interior of the rectangle is the interior of the e-parallelogram. The complement of 
the e-interior in the rectangle is called the e-boundary of the rectangle. See Figure [5721 For 
so € (0,1), the sides of the parallelogram intersect the straight line s = So in two points. 
We denote the ordinates of these two points by g+(so) and g~(so). 
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D 




D 



Figure 4: The e-parallelogram of a rectangle. 



We first state the following useful fact. 

Fact 1. For all 5 G (0, 1), for all A > 0, there exists e, such that for all rectangular boxes 
of side lengths a, b satisfying p = G (5,1 — 6), and all c G [—A, A], the limit shape L P)C 
is entirely included in the e-interior of the box. 

The following lemma controls the function F^} in an e-boundary of a rectangular box 
with side lengths p and 1 — p. 

Lemma 7. Let 5 G (0, 1) and A > 0. Take e as in Fact d Then there exists a positive 
constant C such that for all p G [5, 1 — 6], and all c G [—A, A], 

V( Sj x)G^, Fj$(8,x)<-Cs(l-8), 

where B p is the (macroscopic) box with perimeter 1 and aspect ratio p. 

Proof. For a fixed s, the function x i— > Fp^}(s,x) is concave, and its maximum is reached 
at x = Lp iC (s), and the point (s,L PjC (s)) is inside the e-interior of the box. Therefore 

F$( 8 ,x) < m^{F^(s,gt(s)),F^(s,g-(s))}. 

In a neighborhood of 0, gf(s) = (1 — e)s, and s h-> Fp^}(s, (1 — e)s) X s because it is 
differentiable and vanishes at s = with a non zero derivative. The same argument applied 
to a neighborhood of s = 1, and for g~ gives the result. The uniform bound on C results 
from a compactness argument. □ 
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Lemma 8 (L.U.F.). Let 5 G (0, 1) and A > 0. Take e as in Fact [7J For i/tis e, there exist 
two constants M and k such that if a + b = 2n and £ (5,1 — 5) and for all c G [— A, A], 
then 



-c/n. | ,k. M eX P I ~1ks{\-s) 



V0<k<n, K~b ( X ^k = 2[nL p , c (-) + y^/n~\) < 



\fn ^2-kks(1 - s) 
as soon as (2k, \ nL p ^ c (^) + y-^/raJ ) is in the e-interior of the box. 



Proof. Use Corollary [2j Notice that h c (x,y) x \/~^ uniformly in c. Therefore 



Hp, c (s,x) 



Vp(i-p)^ 



1 



(s - x)(s + x)y (2p - s + x)(2 - 2p - s - x) ' 



The first term is bounded. The second term is of order in a neighborhood of (0, 0) in 
the e-interior of the box. A simple change of variable (s, x) i— > (1 — s, 2p — 1 — x) exchanges 
the second and the third term, which shows that the third term is of order J — in a 

V 1 — s 

neigborhood of (1, 1 — p) in the e-interior of the box. Therefore, 



Hp, c (s,x) 



1 



as long as (s, x) is in the e-interior of the box. 

In order to bound the exponential term in Corollary[2j we bound from below the absolute 

value of 9 F Q C Jf ,X ^ ■ Using that S" (u) = we get 

d 2 F^]}(s,x) _ 2c c c c c 



Qx% gc(s+x) -y gc(s-x) ^ g c (2p— s+x) ]_ gc(2— 2p— s— x) ^ ' 

For s close to 0, the main contribution comes from the second and the third term, that are 
both negative, and of order s _1 . For s close to 1, the main contribution comes from the 
fourth and the fifth term, that are both negative, and of order (1 — s) _1 . □ 

Lemma 9. Let 5 G (0, 1) and A > 0, and e G (0, 1) given by FaciQ Then there exists a 
constant C > such that: for all n > 1, for all c G [—A, A], for all sequences of boxes (B n ) 
with sides a n , b n such that a n + b n = 2n, G (5,1 — 5), 

VA > 0, Vs G [0,1], P(X S G* [-\,\};(X 2ns ,2ns) G Bft < C^L^ll! 
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Proof. Lemma [8] gives the following bound: 



,r 



¥{X S £ [-X,X];(X 2ns ,2ns) G < £ ^f"^ 

l^ A v^ \/2^^(l - i) 

where the index of the sum y runs on the set -^Z — \fnL PjC (s). 

Comparing the sum on the right-hand side and the corresponding integral, we get 



■2 



P(X S [-A, A]; (X 2ns ,2ns) G B° n ) < M' / 2 " s(1 ~ s) da. 

Jr\(-x,x) ^/2ttks(1 - s) 

This last integral is equal to P(UV] > , A ), where Af is a standard Gaussian variable. 

V «*(!-») 

Conclude by using Markov inequality for the fourth moment. □ 
We can now verify Billingsley's condition for tightness [2J. 

Lemma 10. Let ((a n ,b n )) n be a sequence such that a n + b n = 2n, and 2n = P + °(t^)- 
Then there exists a constant C such that for all n > 0, for all < r < s < t < 1, for all 
X > 0, 

K n t(\Xs - X r \ > A; \X t - X s \ > A) < ° {t ~ 02 . (22) 

Proof. Inequality (|22p is automatically satisfied, as soon as \r — s\ or \t — s\ is less than 1/n, 
or A is greater than y/n{\r — s\ A \ t — s\). We suppose now that none of these conditions are 
satisfied. 

Introduce now the three following e-interiors: B £ , the e-interior of the box B := [0, a n ] x 
[0, b n ], B £ L the e-interior of the box Bl '■= [0,ns — ^^2sn] x [0,sn + ^X2 Sn ] and B e R the 
e-interior of the box Br := [ns — \X2 Sn , a n ] x [sn + ^X2 Sn , b n ]. See Fig. [SJ 

Let us first control the probability that the random interface exits the e-interior of the 
box. We use unimodality of the distribution of X 2 nsy Corollary [5J and Lemma [7] to get: 

ns,X 2ns ) B £ ) < 4ens(l - s)F(X 2ns = ng e ± {s)) 

< M ns{1 ~ s) H p . c (s,g £ ± (s)) eM-nCs(l - s)). 
Vn 



Using the uniform bound of H p c >c obtained in the beginning of the proof of 

Vs(l-s) 

Lemma [HJ we get 

M 



P((2ns, X 2ns ) B e ) < M v / ns(l - s) exp(-nC(s(l - s))) < 



(nCs(l-s)) 2 ' 

Using our assumption that n > ^ s _ r i^ t _ s ^i , we obtain: 

\ JTin , Jr - si 4 A \t - s| 4 Ir - si 2 A It - s| 2 „Ji — H 2 , , 

P((2ns, X 2ns ) B £ ) < M A 4 s2(1 _ g) 2 ^ MJ V " M Sv^' (23) 
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Figure 5: The e- interiors B £ (grey), B^f (blue) and (red). 

s 2 A(l-s)'' ! 



Remark 4. The weaker condition n > tkjtz \i yields 



n(2ns,X 2ns ) ? B £ ) < M 5 ^^ 8)2 . (24) 

Now we consider three cases depending on the values of r and t. 
Case 1: r < Then \X r \ < 2r^/n < |. Thus, 

F(\X S - X r \ > A; \X t - X s \ > A) < V££(\X. - X r \ > A) < P(|X S | > ^) 

< P(|X S | > ^,(ns,X ns ) e B e )+F((ns,X ns )?B £ ). 



t - r\ 2 



But we just saw in (J23J) that 

P((2ns,X 2ns ) ?B £ ) < M- 



A 4 ' 

Moreover, according to Lemma [9J 

P(|X S | > ^(2ns,X 2ns ) e B £ ) < C { xi S) < 
Since r < < |, we get that s < |(s — r), and thus 

P(\x s \ > ^,(2n,,x 2ns ) e i? £ ) < < ^j^. 
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And the tightness condition (|22|) is satisfied in that case. 

Case 2: 1 - t < This case is treated similarly as the previous one, and (J22]) is 

again satisfied. 

Case 3: r > -r^= and 1 — t > -A=. This is the generic situation. 

— 4yn — 4yn 

Conditional on {(2ns, A 2ns ) G B £ }, the boxes Bl and Br have an aspect ratio in the 
interval (e, 1 — e). Applying Lemma [9] with 5 = e in the subboxes -B^ and Br, we get 
the existence of an r](e) such that with sufficient probability (2nr, A 2nr ) (resp. (2ni, A 2n f)) 
is in £?2 (resp. B R ^). More precisely, the assumptions on r and t allow to repeat the 
argument to derive (|24|) with the proper scaling in the subboxes i?£ and £>r and get for 
some constant M: 

P((2w, X 2nr ) G* i?f ) |(2ns, A 2ns ) G 5 £ ) < M Ka) < M- — y^~-> (25) 

P((2nt, X 2nt ) ££ (E) |(2n S , X 2ns ) G 5 £ ) < M ^^^f < M^-fil. (26) 

By the Markov property of A, the variables X r and At are independent conditional on 
the value of X s . Hence we can write 

P(|A S -A r | > X;\X t -X s \ > A) = 

= P(A s = fc)P(A r =j|A s = fc)P(A t = /|A s = A:), (27) 

b-fe|>A 
|Z-fc|>A 

where j, k, I runs through all possible values for X r , X s ,Xt (two successive values of j, k, I 
differ by l/^Jn). 

Let E^ t = {(2ns,X 2ns ) E B £ } n {(2nr, A 2nr ) G 5^} n {(2nt, A 2nt ) G From 
Inequalities (|23l) . (|25p . (|26p . there exists a constant M such that: 

ft - r) 2 

i-r(E^ s>t )<M { - 1 p-. 

As a consequence, it suffices to bound the sum in (|27p for values of (j, fc, Z) such that 

{x r = j, x s = k, x t = 1} c {\x s - x r \ > A} n {\x t - x s \ > A} n £* S)t . 

By a scaling argument, when looking at what happens on the left of s, 

Kit & = j\X ns ) = FCj-x^ 2na+x ^ (xr =J- + L p , c (r) - ^TsL ^x 2na J-)) . 

Note that on the right hand side, X r j s is defined with respect to the limit shape inside the 
box [0, 2wa ^ 2 ™ ] x [0, 2rafl +* 2 ™ ]. See Fig. El 
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Figure 6: Limit shapes in the box B and the subboxes Bl and Br. The double arrow 
represents the difference of the limit shapes in B and Bl at time r, whose expression, 
2nL PiC (r) — 2nsL 2ns-x 2ng sc (|), appears in the sticking condition ([28]) . 



When summing over the values of j such that \j — X s \ > A, one can get an upper bound 
for P an ^ n (\X r — X 8 \ > \\X S ) using Lemma[9]as long as the limit shape in the original box 
B and the one in Bl are sufficiently close to each other. We express this proximity by the 
following sticking condition for Bl- 



1 ~ fn 
—j=[X s - A, X s + A] + J-L P)C (r) - y/nsL 2 ns-x 2ns ( 

■\J S y S 2ns ' 



A 



A 



l2s V2s 



(28) 



Indeed, when the sticking condition (128|) is satisfied 
the form: 



, ,lJ n (\X r -X s \ > X\X S ) takes 



■ 2ns — X-2n.s 2ns + ^2 
2 ' 2 



ns ( Ai -7=[^ s - A,X S + A] + \ -L PtC (r) - ^L 2 n S -x 9 , 
s V2s V s 4 ™ s 



-c/n 



< F e 

— r 2"'-^ 2n S + X 2ns 



A A 



2s x/2s 



1 r s— r\2 



< M- 



(\/V2s) 



< 4M 



A 4 



by Lemma [9] in the small box. 

We follow the same idea for the right box Br. Whenever the sticking condition for Br 
is satisfied, 

\t-r) 2 



-c/n 



\X t -X s \>X\X s )<AM- 



A 4 
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Therefore, conditional on one of the sticking conditions to be satisfied, the probability 
of {\X S - X r \ > A} n {\X t - X s \ > A} is bounded by 4M^£. 

Let us now investigate the probability that both of the sticking conditions fail. It is 
bounded by the probability that the one for Bl fails. 

Using Equation (fT4"l) of Lemma[U we replace ^L PtC (r) in (|28l) by L s -L p , c { a ) (§)• 

Condition (|28p is equivalent to 

A 

< -. 

~ 2 



By the mean value theorem applied to p \— > L PiSC (r/s), there exists a p in the interval 
r-f T fr\ T (r\\ 1 ~ dL PjCS (r/s) 



P=P 



But differentiating © with respect to p, for generic values of p, c, i yields 

Q ^ 1 | 1 dL p , c (t) _ e^" 1 ) sinh(c(l - *)) 

2 dp sinh(ct) + e^" 1 ) sinh(c(l - t)) 

„ !0r , i\ sinhfcfl — t)) 1 — i 

sinh(ct) t 

with a constant if which works for all c E [— j4, j4], all p G (0, 1), and all i S [0, 1]. 
As a consequence, the sticking condition (|28p is satisfied as soon as 



\X S \(1 A K --)<-. 



Therefore the probability that (J28J) is not verified is less than 

P(\X a \(l Aif— ) > ^,(2sn,X 2sn ) G 
r 2 

which, by Lemma [9] is bounded by A _4 s 2 (l — s) 2 (l A K A {^ =:L ) i ). This bound is less than 
some constant times A _4 (s — r) 2 < A~ 4 (i — r) 2 , as one can check in both regimes ( — -) *s K. 
This finishes the proof of Lemma [TUJ □ 

6 The unbounded case 

In this section, we relax the constraint of remaining in a box. For any q G (0, 1), we define 
a probability measure ¥ q on all partitions, by 

P«(A) = — i-gl A l,for all partitions A, 
z yQ) 
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where Z(q) = IJfc:i(l ~~ I 1 ) ls the generating function of all partitions. 

When q goes to 1, a law of large number for the shape of the random partition occurs. 
Namely, when rescaled by (1 — q), the boundary of A converges to a deterministic curve. 

Theorem 3 (Vershik). Let L OQ (s) = log(2 cosh s), and (X s ) s ^ be the piecewise linear 
function describing the boundary of X in the (s,x)- coordinates. Then, for all e > 0, 

lim P, (sup |(1 - q)X s(1 _ q) -i - Loo(s)| > e) = 0. 
6.1 Link between L ro and L pc 

Petrov recently discussed in [6] the link between the limit shapes obtained above and the 
infinite shape of the unbounded problem. We explain it in this short section. 

The family of limit shapes obtained above has the following property: if we fix sq 6 (0, 1) 
and take the point (sq, L PjC (so)) as the right corner of a new family of bounding boxes, then 
the limit shape of this new problem is simply the restriction of L P)C to the interval [0, so]; 
rescaled by sq (See Lemma [4]). 




Figure 7: The limit shape in a subbox. 

We are therefore led to the idea that there exists some 'inverse limit' of these curves, 
from which they can all be recovered by restriction. This universal curve is precisely Leo 
from Theorem [3l Indeed, for every p £ (0, 1) and c G M*,, one can find two real numbers 
so < s\ such that the restriction of L ro to [so,si] is (up to an affine transformation) L P)C 
(see Figure [8]) : 

VtG[0,l], L pc (t) = Loo(So + t(Sl " go)) " L - (go) . (29) 

S\ — So 

For given values of so < si, the parameter p is given by 

_ 1 _ -koo(si) - £qo(sq) 

P ~2 2(s!-s ) ' 
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so 



Si 



Figure 8: Limit shapes are restrictions of the universal curve. 



and c is found for example by comparing the slopes of the limit shape in the corner and 



0L o 



(si — so) and 



9L Pz , 



ds 



s c(2p-l)_ cos h c 



s=0 



sinh c 



, which is an increasing function 



identifying 
of c. 

The limit case c — > correspond to si —> sq. For negative values of c one can take 
advantage of the symmetry of the model. 

6.2 Fluctuations around L M 

As in the boxed case, we can compute explicitely the two point distribution. The analogue 
of Proposition [2] in the unbounded case is 

Proposition 5. The 2- dimensional marginal of (X2k)keZ is given by 
[ 2k = 2i,X 2l = 2j) = 



Zk+i,oo(q)Zj+i-i-k,i+i-k-j(q)Zj-i,oc(q) 



,(i-k)(k+i)+(i+j-k-i)U-i) 



where Z a>ao (q) = lim^^ Z a>b (q) = IIj=i(l ~ 9 J ) X - 

The asymptotic analysis of this formula as q — > 1 relies on q-Stirling's Formula (Corol- 
lary [T]) and goes as in Section [5.11 

One can then deduce the following theorem: 

Theorem 4. As q goes to \~ , the random function 

y/2 cosh s 



s i y 



{(l-q)X s(1 _ q) -i -Loo(s)} 
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converges weakly in D to the two-sided stationary Ornstein-Uhlenbeck process (Y s ) s ^, which 
is the Gaussian process on M with covariance 

V(s,t)el 2 , E[Y s Y t ] = e -l*-'l. 

When conditioning the process Y s /\^2coshs to have zero integral over R, one obtains 
a new centered Gaussian process whose covariance is equal to 

e-l f -*l 7T ; 
2 cosh s cosh t 6 

where h(t) = ^tt&nht — log(2cosh£). 

Expressed in the original coordinates (by rotating back the picture by 45 degrees), this 
covariance is the one obtained by Pittel [7] who deals with the microcanonical ensemble of 
partitions (with a fixed area |A| = n — > oo). The fluctuations in that case are thus given by 
this Gaussian, but non-Markov process. 



h(s)h(t), 



7 Appendix: The Ornstein-Uhlenbeck Bridge 

In this short appendix, we give a description of the Ornstein-Uhlenbeck bridge. Let (Bt) 
be a Brownian motion on [0,oo). The Ornstein-Uhlenbeck process (Zt) is the Gaussian 
random process defined by 

Z t = B t -c I Z s ds = e~ ct [ e cs dB s . (30) 
Jo Jo 

Using the right-hand side of (|30|) . we can also represent Z t as: 

Zt = e~ ct B rt eCS ds = e~ rf -B e 2ct_! . (31) 

From (I3ip we derive the covariance. For < s < t, 

E [Z s Z t ] = e" c(s+i) ^^ = e- rfSinh c (C5) . (32) 

'Tying down' this process at t = 1 gives the Ornstein-Uhlenbeck (lt)t6[o,i] bridge of length 
1. Set 

*M = ^ < 33 » 

and define 

Y t = Z t -h(t)Z 1 . (34) 
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Observe that, for < t < 1, 



E[Y t Z x ] =E[Z t Zi-h(t)Zl 



_ c sinh(ct) sinh(ci) _ c sinhc 
c sinh c c 

= 0, 

showing that Yj is independent of Z\. Now let / be any test function and consider 

E [f{X t ) | X 1 = 0] = E [f(Y t + h(t)Xi) | Xi = 0] 
= E [/(Y t ) | Xi = 0] 
= E[/(r*)], 

where the last equality is a consequence of the independence just proved. Thus, the process 
(Yf) has the probability law of the Ornstein-Uhlenbeck process conditional on X\ = 0. Its 
covariance is 



E[Y s Y t ] = E %{Z t - h{t)Z 1 )] = E[Y s Z t 
= E[Z s Z t ] - h(s)E[Z 1 Z t ) 
sinh(cs) sinh(c(l — t)) 



E[Z s Z t ] - h(s)E[Z 1 Z t ] = smh [ cs ] (e" c * sinhc - e - c sinh(ct)) 

csmhc) 



c sinh c 



cf. the middle matrix in (1211) . 

Notice that this covariance is also the Green function of a Brownian motion on [0, 1] 
killed an exponential rate c 2 [3]. Let us mention also, that this process was used by C. 
Donati in her solution of Buffon-Synge problem concerning the typical distance between 
the extremities of a string of given length, thrown at random [4]. 
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